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on the scaling properties of star polymers. We give a brief account of how the 
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1. Star exponents in polymer theory 

It is well known that the asymptotic properties of a flexible polymer chain in a 
good solvent are universal in the limit of infinite chain length. Lowering the temper- 
ature in a system of a polymer in a good solvent the so called temperature may 
be reached below which the polymer collapses. At the transition temperature, the 
effective two point attractive and repulsive interactions between different monomers 
compensate each other and as a result the polymer chain may be described by a 
random walk (up to higher order corrections): the mean square distance between the 
chain endpoint {R 2 ) scales with the number of monomers N as (R 2 ) ~ N. Above the 
G temperature the effective interaction between the monomers is repulsive resulting 
in a swelling of the polymer coil which in the asymptotics again obeys a universal 
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scaling law: 

(R 2 ) ~ N 2u , N^oo (1) 

with the universal exponent v depending on the dimension of space d only. The 
number of configurations Z of a polymer chain scales with iV like 

Z ~ e^PT' 1 , N -> oo (2) 

with a non-universal fugacity e M and universal scaling exponent 7. 

In the early 70-ies, following the work of de Gennes [ |l[ , the analogy between the 
asymptotic properties of long polymer chains and the long-distance correlations of 
magnetic systems in the vicinity of the 2nd order phase transition was recognized 
and elaborated in detail. This mapping allows to define the above exponents ([I]), 
(0) as m — > limits of the correlation length critical exponent v and the magnetic 
susceptibility critical exponent 7 of the 0(m) -symmetric model. In this way, the 
powerful tools of the field theoretical renormalization group approach to the de- 
scription of critical phenomena [ H can be applied in polymer theory (see e.g. [ |3|). 
In particular, this allows to calculate the exponents at d — 3 with record accuracy [ 

v[d = 3) = 0.5882 ±0.0011 and j(d = 3) = 1.1596 ± 0.0020. (3) 

In dimension d = 2, the exact values v — 3/4 and 7 = 43/32 are known due to a 
mapping of two-dimensional (2D) polymers to the 2D Coulomb gas problem [ ||. 

Under similar conditions (i.e. in good solvents above the G temperature and 
long polymer chains), star polymers [ |6j also obey universal scaling laws. For the 
homogeneous star polymer, the asymptotic properties are uniquely defined by the 
number of its constituting chains and by the dimension of space [ |7|. For the partition 
function (number of configurations) Zf of a polymer star of / chains, each consisting 
of iV monomers, one finds: 

Z f ~<? Nf N* l t- x ~(R/e) r, t- fr *, N^oo. (4) 

The second part shows the scaling in terms of the size R ~ N u of the isolated 
coil of N monomers on some microscopic length £, omitting the fugacity factor. The 
exponents 7/, r/f, f = 1, 2, 3, . . . constitute families of 'star exponents', which depend 
on the number of arms / in a nontrivial way. The case of linear polymer chains is 
included in this family with the exponent 7 = 71 = 72 defined in (|2|). 

Apart from rjf and 7/ other sets of star exponents my be defined that govern the 
scaling of different physical properties of polymer stars. However, all these sets of 
exponents can be expressed in terms of a given set by familiar scaling relations [ [F| . 
Moreover, the scaling properties of polymer networks of arbitrary but fixed topology 
are uniquely defined by its constituting stars [|7[], as long as the statistical ensemble 
respects some conditions on the chain length distribution [ |8|] . Thus, the knowledge 
of the set of star exponents 7/ or r]f allows to obtain the power laws corresponding 
to (J|) also for any polymer network of arbitrary topology [ |7|, §|, 0. Recently, the 



118 



Scaling exponents of star polymers 



theory has also been generalized to the case of star polymers and networks that 
contain chains of different species, introducing additional independent families of 
star exponents [ || [11], [1^ . 



As it is usual in the theory of critical phenomena, the space dimension d is crucial 
for the scaling laws (U). At d — 2, there exist exact formulas expressing values of 
the scaling exponents of the star polymers as functions of the number of arms [ 



|10| , 11 , 12] . In d — 3 dimensions on the other hand the problem does not allow for 
an exact solution and different approximate schemes are used. The so-called cone 
approximation takes into account the fact that for large / each chain of the star is 
restricted approximately to a cone of space angle Qf = 47r//. In this approximation 
one finds for large / [ [13|] 

7/ 1- f 2 - (5) 

However, the most accurate values of critical exponents for d = 3 are obtained 
by the field theoretical renormalization group approach refined by a resummation of 
asymptotic series. Just in this way high precision results for the exponents governing 
different critical phenomena are derived [ H , and this is the way the scaling exponents 
(|D for a polymer chain in d = 3 were obtained. In this review, we give a brief account 
on how the numerical values of the star exponents were obtained for d = 3 by means 
of these techniques as well as we provide some examples for phenomena governed 
by the star exponents. 

The paper is arranged as follows. In the next section |2] we give a field-theoretical 
description of a star polymer in d dimensions, the exponents governing scaling of 
a homogeneous star polymer are evaluated numerically in section |3|. The theory 
is generalized to the case of a polymer star consisting of two species of polymers 
(a copolymer star) in section [|. The following sections of the review are devoted 
to the numerical description of different phenomena where star exponents come 
to play: the interaction between polymer stars in a good solvent (section |^), the 
Brownian motion near absorbing polymers (section |6|), diffusion-controlled reactions 
in presence of polymers (section [?[). We finally give some conclusions and an outlook 
in section |8|. 



2. Renormalization of polymer stars 

To give a field-theoretical description of a star polymer let us introduce an Ed- 
wards model of continuous polymer chain [ [3| generalized to describe a set of / 
polymer chains of varying composition, possibly tied together at their end points. 
The configuration of one polymer is given by a path r a (s) in <i-dimensional space lR d 
parametrized by a surface variable ^ s ^ S a . Let us allow for a symmetric matrix 
of excluded volume interactions u a b between chains a,b = 1, . . . , /. The Hamiltonian 
TC is then given by 
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with densities p a (r) — J Q a dsS d (r — r a (s)). In this formalism the partition function 
is calculated clS 8b functional integral: 

Z f {S a } = J P[r a ( S )]exp j-^ft(r*)}. (7) 

Here, the symbol P[r a (s)] includes normalization such that Z{S a } = 1 for all u a b = 
0. To make the exponential of (^-functions in ([?]) and the functional integral well- 
defined in the bare theory, a cutoff so is introduced such that all simultaneous 
integrals of any variables s and s' on the same chain are cut off by \s — s'\ > so- 
Let us note that the continuous chain model (0) may be understood as a limit of 
discrete self-avoiding walks, when the length of each step is decreasing £ — > while 
the number of steps N a is increasing keeping the 'Gaussian surface' S a = N a £^ fixed. 
The continuous chain model (^) can be mapped onto a corresponding field theory 
by a Laplace transform in the Gaussian surface variables S a to conjugate chemical 



potentials ( "mass variables" ) fi a \ |14 



POO 

2fM= / U dS ^ bSbZ d s ^- ( g ) 

Jo b 

The Laplace-transformed partition function Zf{p a } can be expressed as the m = 
limit of the functional integral over vector fields a , a = 1, . . . , / with m components 



a= 1, ...,m : 



Zfi^b} = J V[(f> a (r)) exp[-C{(p b , fi b }] \ m=0 . (9) 

The Landau-Ginzburg- Wilson Lagrangian £ of / interacting fields <pb each with m 
components reads 

W6,M=^E / d ^M' + W a (r)) 2 )+i E u a , a , f d d rcf>l(r)<f>l,(rl (10) 

«=1 J ' a,a'=l J 

here §\ = J^^ii't'a) 2 ■ The limit m — in (0) can be understood as a selection 
rule for the diagrams contributing to the perturbation theory expansions which can 
be easily checked diagrammatically. A formal proof of (|9]) using the Stratonovich- 
Hubbard transformation to linearize terms in (§) is given for the multi-component 
case in [ [TJJ]. The one particle irreducible (1PI) vertex functions T^ L \qi) of this 
theory are defined by: 

S (Z>) r itaM = j eW!d ^i ■ --dr^aM) ■ ••<Mr L ))f PIjW=0 . (11) 

The average (• • •) in ([Ll|) is understood with respect to the Lagrangian ( |I0"D keeping 
only those contributions which correspond to one-particle-irreducible graphs and 
which have non- vanishing tensor factors in the limit m — 0. The partition function 
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Z*f{S a } of a polymer star consisting of / polymers of different species 1, . . . , / 
constrained to have a common end point is obtained from (|7|) by introducing an 
appropriate product of ^-functions ensuring the "star-like" structure. It reads: 

Z> f {S a } = J V[r a ) exp j- J^W(r )} ]j5 d (f a (0) - 5(0)). (12) 

The vertex part of its Laplace transformation may be defined by: 

<J(p + J>)r w) (p,<zi •••?/) = 

e^^d^n . . . d d r f (Mro) ■ ■ ■ M^Mn) ■ ■ ■ 0/(r/))f P T m=0 , (13) 



where all ai,...aj are distinct. The vertex function is thus defined by an 

insertion of the composite operator Y[ a 4>a- Its scaling properties are governed by 
the scaling dimension of this operator. 

In the following we will be mainly interested in two different cases: (i) a polymer 
star consisting of chains of one species with the single interaction u between them 
(i.e. a homogeneous polymer star) and (ii) a polymer star constituted by two species 
of polymers, with interactions Uu, U22 between the polymers of the same species 
and u±2 = 1*21 between the polymers of different species (i.e. a copolymer star). In 
the first case (i) one can also define T^*^ by the insertion of a composite operator of 
traceless symmetry [ [L5[]. In the case (ii) the composite operator in ( |i~3"D reduces to 



the product of two power-of-field operators with appropriate symmetry (0)^(0')^ 2 
each corresponding to a product of fields of the same 'species'. Nevertheless, the 
results are easily generalized to the case of any number of polymer species. 

As is well known, ultraviolet divergences occur when the vertex functions flllD , 
( |T3"D are evaluated naively [ [|. We apply renormalization group (RG) theory to 
make use of the scaling symmetry of the system in the asymptotic limit to extract 
the universal content and at the same time remove divergences which occur for the 
evaluation of the bare functions in this limit [ 0. The theory given in terms of the 
initial bare variables is mapped to a renormalized theory. This is achieved by a con- 
trolled rearrangement of the series for the vertex functions. Several asymptotically 
equivalent procedures serve to this purpose. The results reviewed in the next sections 
were obtained by two somewhat complementary approaches: zero mass renormaliza- 



tion (see [ |2j for instance) with successive e-expansion [ [18| and the fixed dimension 
massive RG approach [ |j~9| . The first approach is performed directly for the critical 
point. Results for critical exponents at physically interesting dimensions d = 2 and 
d = 3 are calculated in an e = 4 — d expansion. The second approach renormalizes off 
the critical limit but calculates the critical exponents directly in space dimensions 
d = 2, d = 3. 

Let us formulate the relations for a renormalized theory in terms of the corre- 
sponding renormalization conditions. Though they are different in principle for the 
two procedures, we may formulate them simultaneously using the same expressions. 
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Note that the polymer limit of zero component fields leads to essential simplifica- 
tion. Each field <p a , mass m a and coupling u aa renormalizes as if the other fields were 
absent. First we introduce renormalized couplings g a b by: 

u aa = fx e Z^Z aa g aa , a = 1,2, (14) 
u 12 = fi £ Z~\ Z~l Z 12 g 12 . (15) 

Here, fi is a scale parameter equal to the renormalized mass at which the mas- 
sive scheme is evaluated and sets the scale of the external momenta in the massless 
scheme. The renormalization factors Z$ a , Z a b are defined as power series in the renor- 



malized couplings which fulfil the following RG conditions: 

d 

Z^igaaj-^^aa^aaigaa)) = 1, (16) 

z l a (9aa)T^ aa (u aa (g aa )) = fi £ g aa , (17) 
Z^ 1 (g ll )Z ( p 2 (g 22 )T i ^ 22 (u ab (g ab )) = ^ £ g l2 . (18) 



These formulas are applied perturbatively while the corresponding loop integrals 
are evaluated for zero external momenta in the massive approach and for external 
momenta at the scale of \i in the massless approach [ [2|, |J . In the massive case the 
RG condition for the vertex function I^ 2 -* reads 

Z^ a (g a a)Taa( U aa(9aa))\k^=0 = /A 1,2. (19) 

In the case of massless renormalization, the corresponding condition reads [ |J : 

ZU9aa)Tf}(u aa (g aa ))\ k 2 =0 = 0, a = 1,2. (20) 
In order to renormalize the star vertex functions we introduce renormalization factors 

z *h,h b y 

ZifzlfZ^ h Y^\u ah {g ab )) = (21) 

With the same formalism one can also describe a star of / mutually avoiding 
walks [ [16[ . In this case all interactions on the same chain u aa vanish and only those 
u ab with a b remain. Then we define the appropriate renormalization for the vertex 
function of mutually avoiding walks (MAW): 

Z f / i 2 Z (MAWf) T* } { AW (u 12 (g ab )) =//'. (22) 

The powers of fi absorb the engineering dimensions of the bare vertex functions. 
These are given by 5f = f(e/2 — 1) + 4 — e. The renormalized couplings g ab defined 
by the relations ([T4]), (|15D depend on the scale parameter /i. By their dependence 
on g ab also the renormalization Z-factors implicitly depend on fi. This dependence 
is expressed by the RG functions defined by the following relations: 

^-T-9ab = Pab(ga'b'), (23) 
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V<t> a (9aa), (24) 

V*fih{9ab), (25) 
vf AW (9a b )- (26) 

The function 77^ describes the pair correlation critical exponent, while the functions 
77*^/2 and rjy AW (g ab ) define the set of exponents for copolymer stars and the stars of 
mutually avoiding walks. Note that Z* 2 o renormalizes the vertex function with a </> 2 
insertion which coincides with r^* 20 \ Consequently, the usually defined correlation 
length critical exponent v is expressed in terms of functions 77*20 and r\<$,. 

Expressions for the (3 and 77 functions will be discussed in the next sections 
together with a study of the RG flow and the fixed points of the theory. 



/i-ln^ /l/2 
//—In Z M AWf 



3. Homogeneous star polymer in three dimensions 

First we give numerical results for a homogeneous polymer star constituted by / 
polymer chains of the same species. The scaling of a star is governed by the laws 
For one species of polymers (a, 6=1), corresponding RG function fl23|)-(|25D depend 
on a single coupling g ab = g and reduce to a single functions (3 a b(ga'b') = (3(g), 
Vj>a(gaa) = V<p(g) and a set of functions r)* hh (g ab ) = r)*(f 1+ f 2 )(g) with f 1 + f 2 = f. 
Here, the function (3(g) is a standard /5-function of the 0(m = 0)-symmetric theory. 
In the infrared-stable fixed point g*\ 

the function rj^g) gives the pair correlation function critical exponent 77: 

V = r](g% (28) 
whereas the functions 77*/) (p) give the scaling exponents 77^ (f|): 

Vf = V*f(9*)- (29) 

The RG functions (f27j)-(|29|) are obtained in the form of perturbation theory series. 
In the massive RG scheme, these are the series in a renormalized coupling g for fixed 
space dimension d whereas the zero-mass renormalization expansion is performed 
both in g and in e = 4 — d. Increasing the order of perturbation theory corresponds 
to increasing the number of loops when the expansion is written in terms of Feyn- 
man diagrams. That is why sometimes the expansions are classified in successive 
approximations in number of loops. The functions (0), (|28| ) are the standard RG 
functions of the 0(m = 0) symmetric theory. They are known by now in high orders 
of perturbation theory, whereas the functions ( |2"§D have been obtained only up to 
three loops. In the field theory they are related to the anomalous dimension of com- 
posite operators of traceless symmetry [ |T5|. In polymer language the star exponents 
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Table 1. Values of the star exponent 7/ obtained in three-dimensional theory 
(columns 2,3,4,5) in comparison with the results of e-expansion (columns 6,7,8) 
and Monte-Carlo simulations (column 9). See the text for a full description. 



/ 


pseudo-e expansion [ 




^-expansion 


MC 




Pade-Borel 


conf. mapping 




m 




mm 


3 


1.06 


1.05 


1.06 


1.06 


1.05 


1.05 


1.07 


1.09 


4 


0.86 


0.86 


.86 


.83 


0.84 


0.83 


0.85 


0.88 


5 


0.61 


0.61 


.58 


.56 


0.53 


0.52 


0.55 


0.57 


6 


0.32 


0.32 


.24 


.22 


0.14 


0.18 




0.16 


7 


-0.02 


-0.01 


-.17 


-.17 


-0.33 


-0.20 






8 


-0.40 


-0.36 


-.63 


-.62 


-0.88 


-0.60 




(-0.99, -0.30) 


9 


-0.80 


-0.72 


-1.14 


-1.11 


-1.51 


-1.01 







were originally calculated to order e in [ |2(|, the e 2 results have been obtained in [ 
0, ^TJ, whereas e 3 expansions are given in [ §. The e-expansion for the exponent rjf 
reads 



Vf = —JU-r 



r f 2 ^77 

l-^(8/-25) + ^[(28/-89)C(3)+8/ 2 -49/+^] 



0{e 3 ) . (30) 



In the massive RG scheme, the three-loop expansions for the d = 3 star expo- 
nents were obtained in the form of pseudo-e expansion in [ F2IJ. However, as already 
noted above, the perturbation theory expansions of the renormalized field theory 
have zero radius of convergence and are asymptotic at best. Therefore, in order to 
make numerical estimates for the expansions one should apply an appropriate re- 
summation procedure. Here, we will not describe the different ways of resummation 
which were applied to the series for the star exponents referring rather to the papers 
I where they are given in detail. For the purpose of the present review we 



just mention that a Borel transform was applied to get rid of the (conjectured) fac- 
torial divergence of the perturbation theory terms [|T7|. Subsequently, the expansion 
for the Borel transformed function was refined by a conformal mapping from the ini- 
tial variable denned on the cut-plane to a disk keeping the origin invariant. Another 
scheme of resummation exploits the analytic continuation of the series for the Borel 
transformed function by means of Pade approximant (Pade-Borel resummation). 

In table [3] we give the numerical results for the star exponents 7/ as obtained by 
the fixed d = 3 RG technique and in the e-expansion in comparison with the data 
of the numerical simulation. First, we give the values of the exponent 7/ obtained 
in the fixed d = 3 technique by Pade-Borel resummation [ Tl\ : the second column 
contains the value of the 7/ obtained directly from the resummation of the series for 
7/, while the third column gives the 7/ based on the resummed series for the star 
exponents x^. These are connected with the 7/ via the scaling relation: 



7/ = l-^/ + K2-W2)-l]/ 



(31) 
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substituting the well-known values of the exponents v{d = 3) = 0.588, r](d = 3) = 
0.027. Columns 4 and 5 give 7/ obtained by the resummation using the conformal 



mapping technique [22]: the resummation of the series 7/(t) (fourth column) and 
the resummation of the series xJt) (fifth column). The next columns show the 
results obtained by the e 3 -expansion using: simple Pade approximation (the 6th 
column) and Pade-Borel analysis neglecting or exploiting the exact results for d = 2 
(the 7th and 8th columns, respectively) [ §]. The last column contains Monte-Carlo 
data [ p5, 26]. For low numbers of arms F ^ 5 the results of the different approaches 



agree reasonably well and are also close to the values obtained by MC simulation. 

The data collected in table [l| was obtained using two different renormalization 
schemes as well as different procedures for the resummation of the resulting asymp- 
totic series. Thus, table [I] gives a test for the stability of the results under the 
changes of the calculational scheme. Obviously, for higher numbers / > 5 of arms, 
the coincidence of the results is no longer good. The main reason for this is that 
calculating the exponents, combinatorial factors lead to an expansion in fe for the 
£ expansion and of fg when directly expanding in a renormalized coupling g. For 
such large values of the expansion parameter, even resummation of the series fails. 
For larger numbers of arms, other approaches to the theory of polymer stars, like a 
self consistent field approximation, might be more useful. 

4. Copolymers and copolymer stars 

Now let us pass to the case of a polymer star constituted by two species of 
polymers. It is described by the Lagrangian (|j~0| ) with interactions Un, U22 between 
the polymers of the same species and Un = U21 between the polymers of different 
species. Such a Lagrangian is used to describe the system of polymers of two species 



immersed in a solvent: the so-called ternary solution [ |27|. A comprehensive analysis 
of the fixed point behaviour of a ternary solution was given in [ [L4| where the RG flow 
of the theory was calculated within the massless renormalization and is known by 
now to the third loop order of the e-expansion. Note, that for the diagonal coupling 
g aa the corresponding expressions are also found in the polymer limit m = of the 
0(m)-symmetric 4 model. They are known in even higher orders of perturbation 
theory. In the massive RG approach, the corresponding expressions were obtained 
in [|J. The equations for the fixed points (FP) P({<?*i, g^i 9n}) °f the /3-functions 
read: 

= 0, a = 1,2, 

/WSri^2*2,<?i* 2 ) = 0. (32) 

As is well known, the first equation has two solutions g* a = Q,g$. For the sec- 
ond equation, one finds a total of 8 FPs depending on the choice for g* a . The 
trivial FPs are G (0, 0, 0), U (# s , 0, 0), U' (0, gfe, 0), S (# s , g£, 0), all corresponding to 
g\ 2 = 0. The non-trivial FPs are found as G(0,0,#q), l%g,0,g{,), U'(0, g$, g^), and 
S(gg, <7g, <7g). In the three dimensional space of couplings #11, #22, #12 these FPs are 
placed at the corners of a cube that is deformed in the gi2 direction (see figure [T]). 
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Looking for the stability of the above 
described fixed points one finds that 
only the fixed point S is stable [14]. In 
the excluded volume limit of infinitely 
long chains, the behaviour of a system 
of two polymer species is thus described 
by the same scaling laws solution of 
only one polymer species. Nevertheless, 
taking into account that real polymer 
chains are not infinitely long, one may 
also find crossover phenomena which are 
governed by the unstable fixed points. 
Knowing the complete RG flow, allows 
one to describe crossover phenomena in 
the whole accessible region [ jl4"| . How- 
ever, for the purpose of our review we 
are interested only in the values of the 
fixed points and in the properties of the 
star vertex functions at these points. 

Let us define the asymptotic values 
of the copolymer star exponents and the 
mutually avoiding walk (MAW) star ex- 
ponents by: 



G U' 



912 



922 



U 



s 



So 



Un 



Figure 1. Fixed points (FPs) of ternary 
polymer solution. The trivial FPs Go, 
Uq, Uq, So correspond to a vanishing mu- 
tual interaction. The non-trivial FPs G, 
U, U', S correspond to a non- vanishing 
mutual interaction (gi2 ^ 0). 



^MAW 



V*flfa(9db)\s, 
V*hf 2 (9ab)\G, 

V*hf 2 (9ab)\u = V*f 2 fi(9ab)\u' 

MAW/ \i 
Vf \9ab)\G- 



(33) 
(34) 
(35) 
(36) 



Taking into account the nature of the fixed points, where the exponents fl3"2"|)-(P6"D 
are defined, one arrives at the following physical interpretation: 

• the exponent describes a star of fi random walks of the 1st type and fi 
random walks of the 2nd type all of which are not selfavoiding; only the chains 
of the 1st type avoid those of the 2nd type and vice versa; 

• the exponent ?7^f 2 describes a star of f\ self and mutually avoiding polymer 
chains walks of the 1st type and f'2 random walks which only avoid those of 
the 1st type; 

• the exponent V^f^ describes a star of fi + f2 random walks all of which are 
not selfavoiding; however, each chain avoids any other chain of the star; 

• exponent V%f 2 describes a star of fi + f 2 self and mutually avoiding chains. 
The exponent can also be expressed by V%f 2 = V^ 1 +f 2 ,o (This relation is broken 
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for the exact results known in d = 2 where a finite effect remains in the zero 
walk limit). This situation has already been discussed in the previous section. 

Starting from the expressions for the fixed points (available both in the form of 
a three-loop e expansion [ [14] or in the form of pseudo-£ expansion at fixed d [ Bl ) 
and the relations ([33|)-(|36l), one can find the series for the star exponents. In the 
e-expansion the following expansions for 77^/3 are obtained [ [|: 

V? lf2 (z) = -/i/ 2 § + /i/ 3 (/ 2 -3 + A)j 

- fx /a (/a - 3 + h) (/1 + h + 3C(3) - 3^, (37) 

vhM = /i(l-/i-3/ 2 )| 

+ ft (25 - 33/i + 8fr 2 - 91/ 2 + 42/i/a + 18/ 2 2 ) ^ 

+ h (577 - 969/x + 456A 2 - 64A 3 - 2463/ 2 + 2290/i/ 2 

- 492/ 1 2 / 2 + 1050/ 2 2 - 504/J2 2 - 108/ 2 3 - 712C(3) + 936/^(3) 

- 224A 2 C(3) + 2652/ 2 C(3) - 1188/ X / 2 C(3) - 540/ 2 2 C(3)) 

(38) 

= -(A-l)/if + /i(/i-l)(2/i-5)^ 

- (A - l)/i(4/a 2 - 20/! + 8/iC(3) - 19C(3) + 25)^. (39) 

Here ((3) ~ 1.202 is the Riemann (^-function. The above formulas reproduce the 
3rd order calculations [ U of the scaling exponents of homogeneous polymer stars 
rjf — rff Q given in formula (^) of the previous section. The corresponding pseudo-£ 
expansion for the exponents 77^/2 obtained in the massive scheme may be found in 
[ |S|. We do not give it here explicitly although in what follows below we will give 
numerical values of the exponents obtained in both approaches. It has been pointed 
out in [ [28| that for the exponent rf^ an exact estimate equal to our first order con- 
tribution may be found. It is indeed remarkable that all higher order contributions 
to rfW vanish in both approaches, rf^ = —£ being an exact result [ ^]. 

With these exponents one can describe the scaling behaviour of polymer stars 
and networks of two components, generalizing the relation for single component 
networks [ 0. In the notation of (|j) one finds for the number of configurations of a 
network Q of Fi and F 2 chains of species 1 and 2 

Z g ~ {Ricr*-***- 1 ***, withr/ g = -dL + J2 N hhVhh, (40) 

h 4-/2^1 

where L is the number of loops and Nf 1 f 2 is the number of vertices with fi and / 2 
arms of species 1 and 2 in the network Q. To receive an appropriate scaling law we 
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Table 2. Values of the copolymer star exponent nf^^ for d = 3 obtained by the 
e-expansion (e) and by fixed dimension technique (3d). 



fl 


1 


2 


3 




4 




5 




6 




h 


e 3d 


e 3d 


e 


3d 


e 


3d 




3d 


e 


3d 


1 


-0.56 -0.58 


-1.00 -1.00 


-1.33 - 


-1.35 


-1.63 - 


-1.69 


-1.88 - 


-1.98 


-2.10 


-2.24 


2 




-1.77 -1.81 


-2.45 - 


-2.53 


-3.01 - 


-3.17 


-3.51 - 


-3.75 


-3.95 


-4.28 


3 






-3.38 - 


-3.57 


-4.21 - 


-4.50 


-4.94 - 


-5.36 


-5.62 


-6.15 


4 










-5.27 - 


-5.71 


-6.24 - 


-6.84 


-7.12 


-7.90 


5 














-7.42 - 


-8.24 


-8.50 


-9.54 


6 


















-9.78 - 


11.07 



Table 3. The values of the copolymer star exponent fn x t 2 at d = 3 obtained by 
the e-expansion (e) and by fixed dimension technique (3d). 



fl 


1 




2 




3 




4 




5 




6 




h 


e 


3d 


e 


3d 


e 


3d 


e 


3d 


e 


3d 


e 


3d 











-0.28 - 


-0.28 


-0.75 - 


-0.76 


-1.36 - 


-1.38 


-2.07 - 


-2.14 


-2.88 


-3.01 


1 


-0.43 


-0.45 


-0.98 - 


-0.98 


-1.64 - 


-1.67 


-2.39 - 


-2.47 


-3.21 - 


-3.38 


-4.11 


-4.40 


2 


-0.79 - 


-0.81 


-1.58 - 


-1.60 


-2.44 - 


-2.52 


-3.33 


-3.50 


-4.28 - 


-4.57 


-5.29 


-5.73 


3 


-1.09 - 


-1.09 


-2.13 


-2.19 


-3.16 - 


-3.30 


-4.20 - 


-4.48 


-5.28 - 


-5.71 


-6.41 


-7.03 


4 


-1.35 - 


-1.37 


-2.61 - 


-2.71 


-3.82 - 


-4.04 


-5.02 - 


-5.40 


-6.24 - 


-6.81 


-7.48 


-8.28 


5 


-1.60 - 


-1.64 


-3.05 


-3.21 


-4.44 - 


-4.75 


-5.80 - 


-6.30 


-7.15 - 


-7.89 


-8.51 


-9.50 


6 


-1.81 - 


-1.89 


-3.46 


-3.68 


-5.01 - 


-5.42 


-6.53 - 


-7.15 


-8.02 - 


-8.92 


-9.50 - 


10.69 



assume the network to be built of chains which for both species will have the same 
coil radius R when isolated. 

In order to calculate numerical values for the exponents ^^j 2 , an d r]y AW in 
[ |[| [29|, |30, [n| the Borel resummation refined by conformal mapping was applied to 
the series fl37D-fl39D as well as to the appropriate series in pseudo-e expansion. The 
resummation showed that the two schemes yield consistent numerical estimates. Our 
tables H, D, and |4] list the results. 

Comparing the numerical values listed in the above tables it is convincing that 
the two approaches and the different resummation procedures all lead to the results 
which lie within a bandwidth of consistency, which is broadening for larger values 
of number of chains /i,/2 > 1. This is not surprising. One can see, e.g. from the 
formulas fl3"7|)-(|3"9]), that the expansion parameters are multiplied by fx and / 2 . 
Rather, it is remarkable that even for a total number of chains of the order of 10 
(see tables |2|, |3|) one still receives the results which are comparable to each other. 

In the next three sections we will briefly review several phenomena where star 
exponents become important. 
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Table 4. The values of the r/^ AW exponents of stars of mutually avoiding walks 
for d = 3, d = 2 obtained by the e-expansion (e) and by the fixed dimension 
technique (3d, 2d). The last column gives the exact conjecture for d 



2 [ 16]. 





d 


= 3 




d = 


2 


/ 


e 


3d 


£ 


2d 


exact 


1 














-.250 


2 


-.56 


-.56 


-1.20 


-1.19 


-1.250 


3 


-1.38 


-1.36 


-2.71 


-2.60 


-2.916(6) 


4 


-2.36 


-2.34 


-4.36 


-4.07 


-5.250 


5 


-3.43 


-3.43 


-6.04 


-5.61 


-8.250 


6 


-4.58 


-4.64 


-7.78 


-7.17 


-11.916(6) 



5. Colloids with polymer stars: The interaction 

Let us now consider the effective interaction between the cores of two star poly- 
mers at a small distance (small on the scale of the size R g of the star). Let us analyse 
the general case of two stars of different functionalities f\ and f%- The power law for 
the partition sum of such two star polymers at a distance r [ |7j, 

zg 2 (r)~r G ^, (41) 

is governed by the contact exponent ©?/ 2 - As we have seen in the section |2] in the 
formalism of the m = component model, the core of a star polymer corresponds to 
a local product of / spin fields 0i(x) • • • </>/(x), each representing the endpoint of one 
polymer chain. The probability of approach of the two cores of the star polymers 
at a small distance r in these terms is described by a short distance expansion for 
the composite fields. The short distance expansion provides the scaling relations 
between the exponents O/^ and the star exponents 7/ (or rjf): 

u ®hh = (7/i - 1) + (7/ 2 - 1) - (7/1+Ja - l ) > 
e /i/2 = Vh + Vfr ~ Vh+h ■ ( 42 ) 

The mean force F?\(r) between two star polymers at a short distance r is easily 
derived from the effective potential Vf^f 2 ( r ) — l°g^/i/ 2 ( r ) as 

Fjl(r) = ^,^ef)~y^. (43) 
The factor 5/18 is found by matching the cone approximation for Q\l (see formula 



(0)) to the known values of the contact exponents for / = 1, 2 [ [34}] . This matching, 
in turn, proposes an approximate value for the 77/ exponents, 

Vf*-^(2 3/2 -2)f/*-f Vl , (44) 



129 



Ch. von Ferber, Yu.Holovatch 



Table 5. The prefactor Qf 1 f 2 of the force between two star polymers at a short 
distance calculated in non-resummed 1-loop and resummed 3-loop RG analysis 
in comparison to the result of the cone approximation. 



©A/ 2 


/2 = 1 


/a = 2 


/ 2 = 3 


/ 2 = 4 


/ 2 = 5 


app. 


/i = l 
/i = l 
/i = l 


-.21 
-.28 
-.27 


-.42 
-.48 
-.45 


-.63 
-.62 
-.60 


-.85 
-.76 
-.73 


-1.06 
-.87 
-.84 


1-loop 
3-loop 
cone 


/i = 2 
/i = 2 
/i = 2 




-.85 
-.82 
-.78 


-1.27 
-1.10 
-1.05 


-1.70 
-1.35 
-1.29 




1-loop 
3-loop 
cone 


/i = 3 
A = 3 
A = 3 






-1.91 
-1.49 
-1.44 






1-loop 
3-loop 
cone 



where —frjx is introduced for consistency with the exact result r\\ = 0. This assump- 
tion nicely reproduces the contact exponents as derived from 3-loop perturbation 
theory [ Q as is displayed in table |5] taken from the [ [33j . 

In table § we use the approximate values of rjf to calculate the cone estimation 
of the contact exponents and compare these with the corresponding values of a 
renormalization group calculation. A perturbation series is taken in terms of an 
e-expansion. The result labelled '1-loop' corresponds to optimal truncation of the 
series by simply inserting e = 1 or d = 3 in the first order term of the expansion. 
The 3-loop result includes a resummation procedure that takes into account the 
asymptotic nature of the series [ [J. The large / result corresponds to the cone 
approximation (|5|). 

The above data convincingly shows that the large / approximation for the short 
distance force between two star polymers can be consistently fitted to the results 
of perturbation theory for low values of the functionality / of the stars. Note, that 
increasing the number of chains /, the polymer stars interpolate between the proper- 
ties of linear polymers and polymeric micelles [ [35| . This approach is general enough 
to describe the interaction between two stars of different functionalities fx and fa. 
This is essential in extending the theory of colloidal solutions of star polymers to 
general polydispersity in / as it appears naturally in any real experiment. 



6. Brownian motion near an absorbing polymer star 

Now, let us turn to another physical phenomenon that is described in terms of 
the star exponents. To this end we consider a model proposed by Cates and Witten 



28 1 to describe the flow of diffusing particles near an absorbing object, which 



may be a polymer chain ( = a selfavoiding walk (SAW)) or a random walk (RW). 
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The flux of particles on the absorbing chain is defined by the density of incoming 
particles close to the absorber (see figure §). Taken p(r) to be the probability density 
of incoming particles, the problem is to solve the steady-state diffusion equation 
(Laplace equation) 



V>(r) = 

with boundary conditions for a given absorbing polymer 



(45) 



p(r) 
p(r) 





Poo 



on the surface of the absorber, 



const 



for 



oo. 



(46) 



One of the motivations for introducing this model was to gain insight for the problem 
of diffusion limited aggregation (DLA) [ [39|. The latter phenomenon is much more 
complicated because of the fact that for DLA the boundary conditions are given on 
the surface of the growing aggregate itself. The process described by equations fl4"5]), 
( fl~6"|) may be rather considered as a diffusion limited catalysis, when particles of one 
type interact with a prescribed fractal (a catalyzing polymer) and transform into the 
other type [ |38fl (see section [7] for more details). 

Let the absorbing polymer of size R 
be chosen from the well defined ensem- 
ble of SAWs. Then one may introduce 
exponents that govern the scaling of the 
moments of the field p(r) close to the 
surface of the absorber (see figure 0). 
For distances R 3> r > a (a being a 
cut off) the averaged moments (p{r) n ) 




are expected to scale as [ ^8 



(p{rT) ~ {R/r)-^\ (47) 

where (...) denotes the average over 
the ensemble of polymers. Taking that 
the flux <p(x) onto any randomly chosen 
point x of the absorber is proportional 
to the field p(r) at a point that is as close as a cut off length a from the absorber, 
one finds for the averaged moments of the flux: 



Figure 2. Star of I random walks at dis- 
tance x from an absorbing polymer of 
size R. 



<0 n > ~ (R/a) 



-A(n) 



(4* 



There is a natural way to associate the exponents A(n) (0), (f4Hj) with the copoly- 
mer star exponents r]f x f 2 034]) , (j35|). The central idea that allows one to calculate the 
properties of the solutions of equation fl4"5| ) with the boundary conditions ([46]) is that 
a path integral representation both for the field p{r) and for the absorber (being a 
RW or a SAW) is possible [ |3SL In terms of the path integral solution of the 
Laplace equation, one finds that p(r) at point r near the absorber is proportional 
to the number of RW that end at point r and avoid the absorber. The nth power of 
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this field is proportional to the nth power of the above mentioned number, i.e. it is 
defined by the partition function of a star polymer with n arms dl2|) (the latter is 
shown by dashed lines in figure 0). Furthermore, introducing the mutual avoidance 
conditions between the "n-arm star" and the 2-arm polymer (representing the ab- 
sorber), one has to calculate the partition function of a co-polymer star consisting of 
chains of two different species that avoid each other. These correspond to the trajec- 
tories of diffusing particles (being RW) and the absorbing polymer (which is chosen 
to be a RW or SAW). Making use of the theory of copolymer stars and networks [ 
[|, one may relate [ the spectrum of exponents (fi7D, ( f£8| ) to the exponents that 
define the scaling properties of co-polymer stars. In particular, as we have seen in 
the section [| for a co-polymer star consisting of f\ chains of species 1 and $2 chains 
of species 2, the scaling of the number of configurations is governed by the expo- 
nents rif 1 f 2 (|40|) . By means of a short-chain expansion [ ^] it is clear that only the 
(smaller) length scale of the absorbing polymer remains and the set of exponents 
Vhf2 can be related to the exponents X(n) (f47[) that govern the scaling of the nth 
moment of the flux onto an absorbing linear chain. Considering the absorber to be 
either a RW or a SAW the exponents read: 

ArwW = ~VL (49) 

AsAwW = -V2U + V20- (50) 



Numerical values for the exponents X(n) may be easily extracted from the nu- 
merical values of the exponents r)2 n given in tables 0. Moreover, one can consider 
an absorber in the form of an /i-arm star and generalize the above considerations 
looking for the scaling laws of diffusing particles in the vicinity of the core of such a 
polymer star. Again, the scaling may be expressed in terms of the co-polymer star 
exponents [ |24|. 



It is well known that polymer chains are fractals with a fractal dimension de- 
fined by a correlation length exponent: df = \ jv. It appears that the behaviour of 
Brownian motion in the vicinity of an absorbing star polymer possesses multifractal 



37fl features. These are analyzed in detail in [ |24l BO 



7. Diffusion-controlled reactions in presence of polymers 

As mentioned in the preceding section, one of the processes described by equa- 
tions (f^), ( f46[) may also be considered as diffusion limited catalysis. In this section, 
we will comment on a general phenomenon which can again be described in terms 
of the star exponents: chemical reactions between diffusing reactants. Examples of 
such reactions can be found in different systems, ranging from biological systems to 



nuclear reactors (see e.g. [ |4(| and references therein). One more place where these 
reactions appear as a limiting stage is the aggregation models [ [|T|. Of particular 
interest are the reactions between reactants of different nature: particles A which 
freely diffuse in a solution and particles B which are attached to polymer chains im- 
mersed in the same solvent [ (the concentration of polymers being low enough to 



132 



Scaling exponents of star polymers 



allow to neglect the inter-chain interaction). Such a process may be also considered 
as a trapping reaction of n particles of A type and traps B: 

A n + B^B. (51) 

The reaction rate k n of ([51]) in the vicinity of a certain trap on the polymer of size R 
is proportional to the averaged moments of the concentration p of diffusing particles 
near this trap. From (f47|) we get that the reaction rate scales with R as [ p8| : 



k n ~ <p»> ~ (R/l)~ Xn , (52) 

with I being a characteristic length scale. 

The scaling exponents A n (|52|) were considered in the previous section (see for- 
mulas (f49|) and (|50|)). Now, considering the absorber to be either a RW star or a 
SAW star of m chains, we define the exponents A mn which can be related to the 
familiar copolymer star exponents r] mn fl35"|) via scaling laws: 



A RW = G 

/y mn Irani 

xSAW _ U , U /co\ 

A mn — 'Iran ^ '/mO' \ °) 

Let us note that the case m = 2 corresponds to a trap located on the chain polymer, 
whereas m = 1 corresponds to a trap attached at the polymer extremity. 

Numerical values for the exponent X mn as well as perturbation theory expansions 
for them follow from tables || ^| and are analyzed in detail in [ ||, |42], |43], |44|. Here, 
let us analyze several particular cases for an absorbing m-star with a reactant trap 
at its core and simply absorbing traps along all chains: 



For a given m-star absorber of size R, the reaction rate ( p2[) scales as k mn ~ 
(R/l)~ Xmn . The increase of the size R by a factor of a results in k' mn ~ 
(aR/iy Xmn leading to: 

k'mn/kmn ~ a~ X ^ (54) 

with X mn positive: increasing R by a factor of a the reaction rate decreases 
a -A m „ tj mes due to the increase of absorbing traps. 

For a given reaction type ([H]) (i.e. for a fixed number n of particles which 
are trapped simultaneously) attaching mi additional arms to an m-arm star 
absorber results in a decrease of the reaction rate: 

hm+m^n/kmn ~ {R/l)-^+^~ X ™\ (55) 

as far as A m2 „ > A min for m 2 > m\. 



For a given m-star absorber the change of the type of reaction ( 51|) to A ni + 
B ^ B, ni > n results in a decrease of the reaction rate: 

kmAn+n^/kmn ~ (iZ/O"^'^)-^, (56) 

as far as A mri2 > A mni for n 2 > m; now, more particles need to find the trap 
at the core simultaneously for a reaction to take place. 
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8. Conclusions and outlook 

The notion of critical exponents is one of the central notions in the theory of crit- 
ical phenomena. The precise determination of the set of critical exponents governing 
an equilibrium thermodynamic second order phase transitions provides a challenge 
in numerous studies. The exponents usually govern the power-law temperature be- 
haviour which manifests itself in the vicinity of the critical point. In polymer theory, 
the direct "relatives" of the exponents governing temperature behaviour at the 2nd 
order phase transitions are the exponents describing the configurational properties 
of long flexible polymer chains. They govern the scaling laws appearing in the limit 
of infinite chain length. These exponents are important in defining the entropy of a 
polymer solution, the radius of gyration of a polymer coil, the osmotic pressure of the 
polymer component in a solution, its rheologic properties, etc [ [I], 0] . Both the 2nd 
order phase transition critical exponents and the polymer chain scaling exponents 
are closely related to each other and, finally, correspond to the second derivative of 
the appropriate free energy. However, in polymer theory there exist sets of exponents 
which do not have direct "relatives" in the theory of 2nd order phase transitions. 
These exponents describe the scaling of star polymers and polymer networks [ [7]. 
Whereas the scaling exponents for chain polymers are calculated by now in different 
techniques with the standard accuracy of modern theory of critical phenomena, this 
is not the case for the star exponents. 

In the present paper we have reviewed some recent calculations of the critical 
exponents governing the scaling behaviour of star polymers. To calculate the criti- 
cal exponents it is today standard to rely on the renormalization group technique 
improved by a resummation of asymptotic series [ |2|, |3|. The numbers provided 
in tables [L]-|4] of the paper provide reliable data which may be used in different 
problems where the star exponents become important. Some of the problems were 
mentioned in the sections [5|-[7]: the interaction between star polymers in a good sol- 
vent, the Brownian motion near absorbing polymers, diffusion-controlled reactions 
in the presence of polymers. 

References 

1. see e.g.: de Gennes P.-G. Scaling Concepts in Polymer Physics. Ithaca and London, 
Cornell University Press, 1979. 

2. Zinn- Justin J. Quantum Field Theory and Critical Phenomena. Oxford University 
Press, 1996; Kleinert H., Schulte-Frohlinde V. Critical Properties of (/> 4 -Theories. Sin- 
gapore, World Scientific, 2001. 

3. des Cloizeaux J., Jannink G. Polymers in Solution. Oxford, Clarendon Press, 1990; 
Schafer L. Excluded Volume Effects in Polymer Solutions. Berlin, Springer, 1999. 

4. Guida R., Zinn Justin J. // J. Phys. A, 1998, vol. 31, p. 8104. 

5. Nienhuis B. // Phys. Rev. Lett., 1982, vol. 49, p. 1062. 

6. We neither give here examples of star polymers nor emphasize their importance in 
applications: rather we refer to the bulk of papers constituting this special issue of 
Condens. Matter Phys. and devoted to the subject. 



134 



Scaling exponents of star polymers 



7. Duplantier B. // J. Stat. Phys., 1989, vol. 54, p. 581-680. 

8. Schafer L., von Ferber C, Lehr U., Duplantier B. // Nuc. Phys., 1992, vol. 374, 
p. 473-495. 

9. von Ferber C, Holovatch Yu. // Phys. Rev. E, 1997, vol. 56, p. 6370-6386. 

10. Duplantier B. // Phys. Rev. Lett., 1986, vol. 57, p. 941-945. 

11. Duplantier B. // Phys. Rev. Lett., 1999, vol. 82, p. 880; Physica A, 1999, vol. 263, 
p. 452. 

12. von Ferber C, Holovatch Yu. // Phys. Rev. E, 2002 (in print). 

13. Ohno K., Binder K. // J. Phys. France, 1988, vol. 49, p. 1329. 

14. Schafer L., Lehr U., Kappeler C. // J. Phys. Paris I, 1991, vol. 1, p. 211-233. 

15. Wallace D.J., Zia R.K.P. // J. Phys. C, 1975, vol. 8, p. 839. 

16. Duplantier B. // Commun. Math. Phys., 1988, vol. 117, p. 279; Duplantier B., 
Kwon K.-H. // Phys. Rev. Lett., 1988, vol. 61, p. 2514-2517. 

17. Hardy G. Divergent Series. Oxford University Press, 1948. 

18. Wilson K.G., Fisher M.E. // Phys. Rev. Lett., 1972, vol. 28, p. 240. 

19. Parisi G. // 1973 (unpublished); J. Stat. Phys., 1980, vol. 23, p. 49. 

20. Miyake A., Freed K.F. // Macromolecules, 1983, vol. 16, p. 1228-1241. 

21. Ohno K. // Phys. Rev. A, 1989, vol. 40, p. 1524-1538. 

22. von Ferber O, Holovatch Yu. // Condens. Matter Phys., 1995, No. 5, p. 8-22; 
Theor. Math. Phys. (Moscow), 1996, vol. 109, p. 34-50 (Engl, ed.: p. 1274-1286). 

23. Semenov A.N. //J. Phys. France, 1988, vol. 49, p. 1353. 

24. von Ferber O, Holovatch Yu. // Phys. Rev E, 1999, vol. 59, p. 6914-6923. 

25. Batoulis J., Kremer K. // Macromolecules, 1989, vol. 22, p. 4277. 

26. Barret A.J., Tremain D.L. // Macromolecules, 1987, vol. 20, p. 1687. 

27. Joanny J.-F., Leibler L., Ball R. // J. Chem. Phys., 1984, vol. 81, p. 4640; Douglas J.F., 
Freed K.F. // J. Chem. Phys., 1987, vol. 86, p. 4280; Schafer L., Kappeler Ch. // 
J. Phys. Paris, 1985, vol. 46, p. 1853; Colloid Polym. Sci., 1990, vol. 268, p. 995. 

28. Cates M.E., Witten T.A. // Phys. Rev. Lett., 1986, vol. 56, p. 2497; Phys. Rev. A, 
1987, vol. 35, p. 1809-1824. 

29. von Ferber O, Holovatch Yu. // Condens. Matter Phys., 1997, No. 10, p. 9-40. 

30. von Ferber O, Holovatch Yu. // Europhys. Lett., 1997, vol. 39, p. 31-36. 

31. von Ferber O, Holovatch Yu. // Physica A, 1998, vol. 249, p. 327-331. 

32. von Ferber O, Holovatch Yu., Schafer L. // Condens. Matter Phys., 1996, No. 7, 
p. 15-25. 

33. von Ferber O, Holovatch Yu., Jusufi A., Likos C.N., Lowen H., Watzlawek M. // 
J. Mol. Phys., 2001, vol. 93, p. 155-158. 

34. Likos C.N., Lowen H., Watzlawek M., Abbas B., Jucknischke O., Allgaier J., Richter D. 
// Phys. Rev. Lett., 1998, vol. 80, p. 4450; Watzlawek M., Likos C.N., Lowen H. 
// Phys. Rev. Lett., 1999, vol. 82, p. 5289; Jusufi A., Watzlawek M., Lowen H. // 
Macromolecules, 1999, vol. 32, p. 4470. 

35. Gast A.P. // Langmuir, 1996, vol. 12, p. 4060; Grest G.S., Fetters L.J., Huang J.S., 
Richter D. // Adv. Chem. Phys., 1996, vol. 94, p. 67. 

36. von Ferber C. // Nucl. Phys. B, 1997, vol. 490, p. 511-542. 

37. Hentschel H.G.E., Procaccia I. // Physica D, 1983, vol. 8, p. 435-444; Halsey T.C., 
Jensen M.H., Kadanoff L.P., Procaccia I., Shraiman B.I. // Phys. Rev. A, 1986, 
vol. 33(2), p. 1141-1151. 

38. Burlatsky S.F., Oshanin G.S., Likhachev V.N. // Sov. J. Chem. Phys., 1991, vol. 7, 



135 



Ch. von Ferber, Yu.Holovatch 



p. 1680-1698; Burlatsky S.F., Oshanin G.S. // Phys. Lett. A, 1990, vol. 145, p. 61-65; 
Oshanin G., Moreau M., Burlatsky S. // Anv. Coll. Int. Sci., 1994, vol. 49, p. 1-46. 

39. Witten Jr., T.A., Sander L.M. // Phys. Rev. Lett., 1981, vol. 47, p. 1400-1403; Phys. 
Rev. A, 1983, vol. 27, p. 5686-5697. 

40. Rice S.A. Diffusion-limited reactions. - In: Bamford C.H., Tipper C.F.H., Comp- 
ton R.G. (eds). Comprehensive Chemical Kinetics. Vol. 25. Amsterdam, Elsevier, 1985. 

41. Meakin P. - In: Domb C, Lebowitz J.L. (eds). Phase Transitions and Critical 
Phenomena. Vol. 12. New York, Academic Press, 1988, p. 336-442; Simulation 
of Non-equilibrium Growth and Aggregation Processes. Landau D.P., Moon K.K., 
Schiittler H.-B. (eds). Springer Proceedings in Physics. Computer Studies in Con- 
densed Matter Physics. Vol. 33. Berlin, Springer- Verlag, 1988. 

42. von Ferber C, Holovatch Yu. - In: "Renormalization Group'96" 3rd Int. Conf., Dubna, 
Russia, August 1996, Ed. by: Shirkov D.V., Priezzhev V.B. JINR, Dubna, 1997, p. 
123-138. 

43. von Ferber C, Holovatch Yu. // Condens. Matter Phys., 1999, vol. 2, p. 303-312 

44. von Ferber C, Holovatch Yu. // J. Mol. Phys., 2001, vol. 93, p. 151-154. 



136 



